Entanglement scaling in two-dimensional gapless systems 
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We numerically determine subleading scaling terms in the ground-state entanglement entropy of 
several two-dimensional (2D) gapless systems, including a Heisenberg model with Neel order, a free 
Dirac fermion in the 7r-flux phase, and the nearest-neighbor resonating-valence-bond wavefunction. 
For these models, we show that the entanglement entropy between cylindrical regions of length x 
and L — X, extending around a torus of length L, depends upon the dimensionless ratio x/L. This 
can be well-approximated on finite-size lattices by a function ln(sin(7ra::/L)), akin to the familiar 
chord-length dependence in one dimension. We provide evidence, however, that the precise form of 
this bulk-dependent contribution is a more general function in the 2D thermodynamic limit. 



Introduction - The study of quantum condensed mat- 
ter systems is benefiting from an infusion of ideas related 
to quantum information and entanglement. The impor- 
tance of this new resource is strikingly demonstrated in 
the study of entanglement entropy at one-dimensional 
(ID) quantum critical points with conformal invariance. 
Conformal field theory (CFT) provides an important uni- 
versal number, the central charge c, which appears in an 
astonishing array of physical quantities.^ A given CFT, 
and thus any quantum critical points it describes, can be 
characterized by this number. Its numerical or analytical 
determination provides an invaluable tool for identifying 
which, if any, CFT describes the scaling limit of a given 
Hamiltonian. Computing the entanglement entropy has 
proven to be a very useful way of finding c numerically. 
It can be extracted directly from the ground-state wave- 
function by measuring its Renyi entanglement entropy, 
Sn = 1/(1 — n)ln [Trp^], where region A is entangled 
with its complement, region B. Namely, in a system 
with total length L, where the region A has length x, 
the scaling of the Renyi entropy in ID critical systems 
depends on the "chord length" as,!^^ 

^n = ^ 1 + - In -sin — , 1 

with the central charge appearing as the coefficient. 

In higher dimensions, the scaling behavior of the entan- 
glement entropy is much less well understood. Ground 
states of local Hamiltonians are generally believed to pro- 
duce an "area-law" (i.e. boundary) scahng,'^] the sublead- 
ing corrections to which may be universal quantities that 
can be used to identify and characterize quantum phases 
and phase transitions. A well-established example of this 
is the topological entanglement entropj^l^ of a gapped 
state with topological order. In gapless states, the sub- 
leading corrections may still potentially harbor universal 
quantities. It is conceivable that such quantities could be 
used to define an "effective" central charge in two spa- 



tial dimensions, but there are strong constraints on any 
proposal.^* The best-understood gapless situation in two 
dimensions is the special case of a conformal quantum 
critical point, where the ground state itself is written in 
terms of a two-dimensional (2D) CFT.E^tl^ In the pres- 
ence of a spontaneously broken continuous symmetry, 
Goldstone mod es produce a subleading bulk logarithmic 
correction.EsHni Subleading logarithms from corner con- 
tributions wit h unive rsal coefficients also occur at some 
critical points.'i^EniS! 

The purpose of this paper is to analyze one type of sub- 
leading term in 2D gapless systems and to study whether 
this term is universal. Gapless modes typically have long- 
range correlations, so it is possible for the entanglement 
entropy to depend on the size and shape of regions A and 
B. Indeed, the ID result, ([T]), is manifestly size depen- 
dent. We show how similar behavior also occurs in two 
deimensions. 

We study the finite-size scaling of the second Renyi 
entropy for the ground states of several 2D gapless sys- 
tems on the square lattice using quantum Monte Carlo 
(QMC) simulations. It is possible to vary the size of 
regions A and B without changing the length of the 
boundary between a toroidal lattice geometry, where A 
and B are cylinders as in Fig. [T] We examine the Neel 
ground state of the Heisenberg model, and the nearest- 
neighbor resonating-valence-bond (RVB) wavefunction, 
in this geometry. In both cases, we find a size- and 
shape-dependent scaling function that closely mimics the 
chord-length contribution in one dimension in Eq. ([T]). 

To probe this behavior in a simpler system, we also 
study free spinless fermions in the 7r-fiux phase and find 
that the entanglement scaling also has a universal size- 
and shape-dependent piece. For finite-size systems, this 
closely mimics the chord length, but in the infinite-size 
limit we observe it to cross over to a different nontrivial 
function. Among other consequences, this term will give 
a nonzero signature in the entanglement quantitiegSHSl 



FIG. 1: (Color online) An 8 x 16 toroidal lattice. The width 
of cylindrical region A (blue) is x = 4. The boundary length 
between region A and its complement is ^ = 16. 



designed to look for topological order, which complicates 
any possible generalization of the topological entangle- 
ment entropy to gapless spin-liquid states. 

Fermions with it flux — We begin by considering free 
spinless fermions on a square lattice, with n flux through 
each plaquette. We consider a torus of size by Ly, 
and measure the entanglement using a cornerless cylin- 
drical region A (Fig. [I]) with a constant boundary length 
£ = 2Ly. We denote the width of region ^ by a;. This 
system has Dirac points near momentum ky = and 
ky = IT. We take antiperiodic boundary conditions in the 
X direction so that there will be no exact zero mode. We 
use exact numerical diagonalization of the single-particle 
Hamiltonian to compute the entropy. The entanglement 
entropy of 2+1-dimensional conformally invariant sys- 
tems such as this has been argued to be of the formj^^^HH 



Sn ~ const. X e/a + ^{x/L^, Ly/L^), 



(2) 



where 7 is a universal scaling function of the dimen- 
sionless ratios. The area-law term proportional to the 
boundary length £ depends on the lattice constant a, 
and so the constant is nonuniversal. A crucial differ- 
ence from the result in one dimension, Eq. ([I]), is that 
a only appears in the area-law term. In contrast, the 
ID result can be written as the sum of two terms, as 
Sn = Cln[sin(f;)] + Cln[^], where C = c/6(l + 1/n). 
The first term is a universal function of the dimensionless 
ratio x/L, akin to the function 7 above, while the second 
term involves the lattice scale, as it diverges with L. 

To illustrate the absence of such an "additive loga- 
rithm" (a logarithmic divergence depending on L/a) in 
two dimensions, we treat this free system as a collec- 
tion of independent systems in one dimension labeled 
by the momenta ky. The ky — mode contributes an 
additive logarithm Cln(L^) to the entropy, while the 
modes wi th sm all ky ^ Q contribute additive logarithms 
C In ( fc^ ^ ) !21ilSI Summing over ky — 211] /Ly, this gives an 

entropy C[ln(L,) + 2EjC['' ln(L;,/27rj)] = Cln(i,) + 



2Cln[(L3,/27r) 



where the factor of 2 arises from 



summing over positive and negative m ^ 0. Using Stir- 
ling's formula for Ly\, one finds that the additive loga- 
rithm terms add to C[ln(L2;) — ln(Lj,)] = Ch\{Lx/ Ly). 
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FIG. 2: (Color online) (a) Renyi entropy of the ~ 256 
Dirac fermion model, with an arbitrary constant subtracted 
off each data set. The straight line is Eq. ([l]) plotted with 
c = 2. It is clear that deviations from linearity increase for 
larger Ly/L^. (b) A set of points near a; = is shown to 
follow the ID result ln(sin(7ra;/I/3;)). Note that the "central 
charge," c, in this region is independent of the system size. 



This can be absorbed into the scaling function 7, so that 
there is no additive logarithm. A more precise calcula- 
tion would include the effect of finite L^, but we ignore 
this since it does not affect the cancellation of additive 
logarithms. A similar calculation near ky — n leads to a 
cancellation of the additive logarithm there. 

The entropy of a given ky mode contains, in addition 
to the additive logarithmic divergence in ky, a universal 
scaling function G{x / Lx,kyx). At ky = 0, we see the 
chord-length scaling Cln [sin(^)] (Fig. [2]), but for ky ^ 
and for kyX large, the chord-length scahng disappears 
and the entropy becomes roughly flat as a function of 
x/L. In fact, for Ly — L^ ~ L, the lowest ky mode 
has a mass 2ky = An/L. This factor of 47r « 13 means 
that this mass is rather large, and so the entropy of this 
mode is flat for a large range of x/L. As a result, for 
Ly = Lx = L, the entropy of the 2D system appears 
to display ID chord length scaling over a wide range of 
x/ L,j-. 

Quantum Monte Carlo — Using QMC techniques, we 
simulate both the Heisenberg ground state and the RVB 
wave function in two dimensions. The Heisenberg ground 
state is projected from a trial state by applying a high 
power of the Hamiltonian, H — ' '^i^ ^ QMC 

method operating in the valence bond (VB) basis. The 
RVB wavefunction is an equal-amplitude superposition 
1^) = 12a 1^") nearest-neighbor VB states, 

N/2 

1^") = o]v74n(lt4,J-|U,J), (3) 



1=1 



defined by requiring that each spin i on one sublattice be 
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FIG. 3: (Color online) The second Renyi entropy for the Neel 
and RVB states for L — 24. Note that the entanglement 
entropy for the RVB splits into two branches, even and odd, 
which may be related to the existence of topological sectors 
in the underlying transition graphs. 
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FIG. 4: (Color online) (a) Heisenberg data and linear fits 
(excluding the first two data points on the left) for L = 
14, 16, . . . , 28 plotted in terms of the log of the "chord length," 
In [sin^]. (b) Slopes of the fits, c(L), exhibit a strong de- 
pendence on the system size, L. 



in a singlet with one of its nearest neighbors Jq, .I^^^^ZI The 
RVB Monte Carlo sampling algorithm does a random 
walk through the possible states by creating a defect at 
some spatial point and propagating it through the system 
(thereby rearranging the nearest-neighbor bonds) until 
the defect reaches the initial point and its path forms 
a closed loop-t^^l If we visuaUze the Heisenberg ground 
state in this VB language, then the RVB wave function 
is its largest component, the remainder of the state be- 
ing superpositions of longer bonds, decaying with their 
length as Xjr^ ?^ Likewise, the RVB wave function is the 
ground state of a local (but longer range) Hamiltonian 
that includes a Heisenberg term.l^Sl 

We consider the same geometry as for the vr flux 
fermions, with = Ly = L (see Fig. [ij. In Fig. |3j 
we plot QMC results for the second Renyi entropy in 
the Neel and RVB states on a 24 x 24 torus. Several 
features of the entanglement scaling are clear from this 
plot. First, note that the data for the Neel state has 
a significant curvature as a function of x. This curva- 
ture was first seen in Ref. |19j but not explored in detail 
[instead, using a fixed cc/L, a surprising subleading log- 
arithmic term oc \a{i) was founcP^. The entropy of the 
RVB wave function exhibits an obvious dependence on 
whether x is even or odd which we discuss in more detail 
below. In each of the even and odd "branches," there 
is significant curvature in the x dependence as with the 
Heisenberg case. 

To capture the x-dependent curvature of these wave 
functions, we fit the data with the scaling ansatz. 



S2 = al + b\n{l)+c{L) In 



fiix\ 



(4) 



motivated by the chord length in Eq. ([T]). We begin by 
examining the Neel state in Fig. |4j For a fixed linear 
system size L and boundary length £, plots of 5*2 versus 



In [sin(^)] would yield a straight line if Eq. (W| were 
obeyed perfectly. The plots indeed are quite close to 
straight lines for a fixed L. 

The second Renyi entropy therefore displays, at the 
very least, an effective chord- length dependence over a 
large range of x for the square torus. It is possible that 
the apparent chord-length scaling of this 2D system is not 
perfectly obeyed in the thermodynamic limit and that 
this fact is manifest in slight deviations from straight-line 
behavior in Fig. Qa) . This would be a similar scenario 
to the deviation from chord-length scaling observed for 
TT-fiux fermions in Fig. [2j However, it is difficult to draw 
a firm conclusion regarding the statistical significance of 
any deviation from Eq. (|4| scaling in our present data, 
due to limited system sizes and stochastic error. 

Wc can, however, further examine the deviation from 
conformal-style scaling by extracting the L dependence 
of the coefficient c{L) in Eq. Q. In order for this shape- 
dependent term to be universal in two dimensions, c{L) 
should approach a constant in the limit L ^ 00 for fixed 
x/L. As illustrated in Fig. |4|^b), the coefficient does not 
approach a constant for the system sizes that we have 
studied but, rather, has some functional dependence on 
L. This functional dependence is apparently sub linear 
- possibly behaving like c(L) ~ with p < 1. This 
scenario could be supported by the QMC data if con- 
vergence were assumed to be very slow. Indeed, in the 
quantum dimer model, the corresponding term can be 
computed exactly in finite size, and the convergence is 
very slow.l^ A definitive determination of this limit (and 
therefore the strict adherence of 7 to universality in this 
system) is impossible with our current data; significantly 
larger system sizes must be studied. 

We next examine the scaling of the Renyi entropy in 
the RVB wave function. As shown in Fig. [3j a strik- 
ing two-branch structure exists, depending on whether 
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FIG. 5: (Color online) (a) Even and (b) odd branches of the 
second Renyi entropy plotted against the log of the "chord 
length." We exclude x = 1,2 data from the plots, as there is 
some crossover behavior shown in Fig. 3] (c) Absolute value 
of the slopes, |c(Z/)|, as a function of the system size, L. As 
with the Heisenberg, there is a strong dependence on L. 



the distance x is even or odd. The presence of the two 
branches presumably is related to the fact that correla- 
tors in the RVB state have a pronounced even-odd depen- 
dence. Moreover, simple counting arguments of prototyp- 
ical VB configurations in the (0,0) topological sectop25EZl 
show that the number of VBs crossing from region A to 
region B alternates strongly with x. This L = 24 data 
display a clear a;-dependent curvature in each branch. 
This can be analyzed more closely by attempting fits of 
the form of Eq. (|4| to each branch individually. From 
Figs, [sja) and Sjfb), it is clear that fits to the scaling 
ansatz to both branches are quite accurate when the ex- 
tremal values of x are excluded. It is worth noting that in 
the closely related quantum dimer model on a square lat- 
tice, similar terms appear and can be computed exactly,'^ 
generalizing the results of Refs. [TSHIH]. 

We can attempt to extract the size dependence of the 
coefficient c in a similar manner as for the Heisenberg 
model. Unfortunately, due to the two-branch structure, 
each curve in this plot has essentially half the usable data 
compared to the analogous Heisenberg results in Fig. [4] 
Nonetheless, the result [Fig.jsjjc)] shows that a significant 
L dependence seems to exist in the RVB wave function 
as well. This again suggests that, although the fit to a 
chord-length scaling at fixed L is consistent within the 
accuracy of our data, subtle corrections to this form may 
come into play in the 2D thermodynamic limit. 

Discussion — We have studied the Renyi entanglement 
entropy in the ground state of three gapless systems on 
Lx X Ly toroidal lattices, where the subregion A is a 
cylinder of length x. We have demonstrated that it 
contains a subleading scaling term which depends on 



bulk quantities, namely the dimensionless aspect ratios 
of the subregion and the lattice linear dimensions. Note 
that while num erica l measurement of topological entan- 
glement entropy'^'^^' has been used to probe topological 
properties of gapped phases,'^ the subleading term con- 
sidered here means that a measurement of topological en- 
tanglement entropy in a gapless phase could give either 
a zero or a nonzero result, even without any topological 
aspects of the phase [though measurements in the C/(l) 
superfluid phase yielded a vanishing numbeil^ . Interest- 
ingly, just as strong subadditivity constrains the sign of 
the Levin- Wen entropy,^ it also implies, for fixed Lx,Ly, 
that 7 for the von Neumann entropy is a concave-down 
function of x. 

Our QMC simulations of the Heisenberg Neel ground 
state and the short-range RVB wave function with — 
Ly ~ L show an almost-perfect logarithmic dependence 
of 7 on the chord length sin(7ra:/L). It appears that the 
coefficient of this term is not a universal constant, how- 
ever, which might suggest either that care must be taken 
in the order of limits with which the thermodynamic limit 
is approached or that a size dependence remains in this 
limit, rendering this term nonuniversal. A study of the 
crossover from one to two dimensions might illuminate 
this issue further. Further evidence that the true 2D 
scaling function might not be exactly the chord-length 
form is given by the scaling of gapless Dirac fermions in 
the vr-flux phase. Here we have argued that such scaling 
is superseded by a sum over transverse modes, leading 
to a different (unknown) functional form in two dimen- 
sions. Furthermore, spontaneous symmetry breaking in 
the Heisenberg model may complicate measurement of 
the entanglement entropy. The fact that a complete char- 
acterization of the scaling behavior in the Neel and RVB 
states remains a challenge, despite the large lattice sizes 
studied to date, underlines the absolute necessity for us- 
ing large-scale QMC simulations for the study of entan- 
glement entropy. 

Regardless of the precise functional form of the shape- 
dependent subleading term 7, its general existence in 
gapless wavefunctions in two dimensions would have 
some profound consequences. Besides the immediate 
complications in attempting to use entanglement as a 
probe to detect gapless spin liquids mentioned above, 
the similarity of the scaling function to a chord length 
(present in ID conformally invariant systems) raises the 
tantalizing possibility that our results will prove useful in 
characterizing higher dimensional critical points. Indeed, 
since the search for a c theorerrP^ valid in higher dimen- 
sions is o f inte nse interest across several disparate field 
of physics, l^^'^^Hnil -^^g hope that our results will inspire a 
broader examination of this scaling term in 2D gapless 
states. 
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